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Abstract 

We devise a hierarchy of computational algorithms to enumerate the microstates of a system 

comprising A^ independent, distinguishable particles. An important challenge is to cope with 

^vq ' integers that increase exponentially with system size, and which very quickly become too large to 

be addressed by the computer. A related problem is that the computational time for the most 
^ ■ 
^~-*', obvious brute-force method scales exponentially with the system size which makes it difficult to 

study the system in the large A^ limit. Our methods address these issues in a systematic and 

hierarchical manner. Our methods are very general and applicable to a wide class of problems 

Y^ . such as harmonic oscillators, free particles, spin J particles, etc. and a range of other models for 

• 1—1 ' 
£^ ' which there are no analytical solutions, for example, a system with single particle energy spectrum 



given by e{p, q) = eo(p^ + Q'^); where p and q are non-negative integers and so on. Working within 
the microcanonical ensemble, our methods enable one to directly monitor the approach to the 



K*" I thermodynamic limit (A^ -^ oo), and in so doing, the equivalence with the canonical ensemble is 

_^ ' made more manifest. Various thermodynamic quantities as a function of N may be computed using 

. ■ our methods; in this paper, we focus on the entropy, the chemical potential and the temperature. 
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I. INTRODUCTION 

Entropy is critical to the understanding of statistical physics. However, it is our expe- 
rience that students have difficulty in conceptualizing entropy. For example, the follow- 
ing simple problem has been presented to introductory statistical physics students at our 
institutions^ on numerous occasions over more than a ten year span, and the question has 
usually been bewildering to many of them: 

Imagine a container with a large number of red marbles and another container with a similar 
number of green marbles. Mix the marbles together, and shake the container really well. After 
this, a handful will draw, on average, an equal number of red and green marbles. Question: 
How does this come about - do the red marbles 'know' of the green marbles, and vice versa?! 



Of course the red marbles do not communicate with the green marbles! The answer has 
to do with the fantastically large number of ways in which the system can configure itself 
for which the final state may be regarded as a mixed state compared to a separated or even 
partially separated state. 

Computation has a significant role to play in enhancing the understanding of entropy. 
Entropy, as in the above example, is often presented in the context of the second law of 
thermodynamics as a measure of the state of randomness of a system. The Monte Carlo 
method, therefore, is a very obvious computational method in the study of entropy, and many 
pedagogical texts focus on the use of the random variable to simulate various thermodynamic 
pro cesses .-"- 

In this paper, however, we approach entropy by computing the relevant statistical quan- 
tities more directly. Our starting point is the enumeration of system microstates. This 
very quickly becomes a computational challenge because we have to cope with integers that 
increase exponentially with the system size, and therefore become larger than the largest 
integer that can be addressed by a computer. A ffist step in our hierarchy of computational 
methods to address this problem is to work with real numbers using real arithmetic, but 
this improvement is minuscule compared with the scale of the problem. A related problem 
is that the computational time for the most obvious brute-force method scales exponentially 
with the system size which makes it difficult to study the system in the large N limit. 



We devise a hierarchy of computational algorithms to enumerate the microstates of a 
system comprising A^ independent, distinguishable particles. Our methods are very general 
and applicable to a wide class of problems such as e{p) = EqP (quantum harmonic oscillators), 
e{p) = Eop"^ (ID free particles), e{p,q,r) = £o(p^ + 0'^ + r^) (3D free particles), spin J 
particles with e{m) = Eqitl where m = — J, — J + 1,--- ,J, etc. and also a range of other 
models for which there are no analytical solution, for example, e{p) = e^p^, e{p, q) = 
eo(p^ + Q^), etc. Here p, q and r are integers. J is integral or half odd integral for real spin 
systems. The systems may be physically realizable or they may be idealized. They may 
or may not be solvable analytically. These models may have a finite or an infinite single 
particle energy spectrum. Our numerical method is applicable to all of these systems and 
has been tested against various analytically solvable models such as a system of harmonic 
oscillators, a system of free particles, a system of spin 1 particles, etc. This has enabled 
us to check the accuracy of our numerical methods. Working within the microcanonical 
ensemble, our methods enable one to directly monitor the approach to the thermodynamic 
limit (A^ — 7- oo), and in so doing, the equivalence with the canonical ensemble is elucidated. 
Various thermodynamic quantities as a function of A^ may be computed using our methods; 
in this paper, we focus on the entropy, the chemical potential and the temperature. 

Furthermore, our work elucidates different computational algorithms that give the student 
an experience of different ways of solving the same problem, with some methods being 
exponentially more efficient than others. This point is worth pondering about: We do 
not seem to be spending sufficient time focusing students' minds on the important but 
distinct aspect related to algorithm development in the teaching of computational physics. 
Very often, we go from discussions of a physics problem directly to code development, often 
skimping on the details of algorithm development. Our presentation in this paper encourages 
the student to think a little more carefully about the manner in which the problem should 
be solved, and gets the student to appreciate how a poorly thought algorithm, even though 
mathematically correct, can lead to catastrophe. 

In our work, temperature as a statistical quantity emerges very naturally. For small 
system sizes, the curve of entropy versus energy is jagged and thus not differentiable, which 
makes the temperature of the system ill-defined since 
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As the large system size limit is approached, we show that the curve of entropy versus 
energy converges to a smooth, differentiable form for which temperature is numerically 
better behaved; this is consistent with the fact that temperature for an isolated system only 
really makes sense in the thermodynamic limit. Our approach gives a clear indication of the 
rate at which the thermodynamic limit is reached, and what 'large' actually means from a 
computational point of view. 

The subject of the enumeration of states is, of course, an established one^ that has been 
developed since the time of Boltzmann, more than a hundred years ago. However, publi- 
cations over the course of the past ten years^"— show that the pedagogical understanding 
of the subject is still in need of attention, and whilst the theoretical concepts are well un- 
derstood, the numerical evaluation of results is still lagging behind arguably for the reasons 
given above, namely, the numbers increase too rapidly which make them difficult to handle 
comput at ionally. 

Analytical expressions for the number of system microstates, namely, Q{N, E) for A^ 
particles with total energy E, only exist for a small number of problems such as a system 
of A^ independent quantum harmonic oscillators or a system of A^ free particles. Here the 
Stirling approximationii is applicable and enables the derivation of closed form analytical 
results. For a wide range of other problems, such analytical solutions do not exist, and 
accurate numerical methods are needed. 

Moore and Schroeder- considered a system of quantum harmonic oscillators and noted 
that the results become 'rather cumbersome with much more than a few hundred energy 
units, while [computer] overflow errors can occur if there are more than a few thousand 
oscillators ...'. We present in this paper a robust means of avoiding this type of overflow 
error. Prentis and Zainiev^ considered systems in the large A^ limit. In our work, we compute 
Q{N,E) recursively for A^ = 1, 2, 3 ■ ■ ■ N^ax where N^ax can be made sufficiently large to 
test various results in the large A^ limit, i.e using our method, we are able to approach 
the thermodynamic limit systematically and with no approximations. Our algorithm is 
numerically stable, and the results are achievable in real-time on a moderate modern desk- 
top computer.— Schoepf^ considered equally spaced energy levels in studying entropy and 
the approach to thermodynamic equilibrium, and remarked that 'students have a better 
intuitive feel for the concept of energy, an intuition that is lacking for entropy', a point that 
we strongly agree with. 



In Section II, we present a non-trivial model for which no analytical results exist. We 
discuss a hierarchy of methods that enable a complete and numerically exact resolution of 
the problem. Working within the microcanonical ensemble, we give results for the entropy 
and we compare our results as a function of A^ with the equivalent results in the canonical 
ensemble which are, by definition, in the thermodynamic limit. This enables us to directly 
monitor the systematic approach to the thermodynamic limit and to test the equivalence 
of the two ensembles. This has enormous pedagogical value. We consider the temperature 
and the chemical potential, and we make similar comparisons with equivalent results in the 
infinite size limit. In Section III, we draw our conclusions and give suggestions for further 
work. 

II. THE COUNTING OF MICROSTATES 

A. The model 

As indicated in the introduction, our methodology is very general and applicable to a 
very wide range of different problems. Here we have chosen to focus on the following model 
because it is not analytically solvable and because it has non-trivial structure. This model 
is more complex than the analytically solvable models. We consider an isolated system 
consisting of A^ identical, non-interacting quantum particles with single-particle energy given 
by 

e{p,q) = eo{p^ + q'), (2) 

where p and q are non-negative integers. We wish to determine the total number of system 
microstates Q{N, E) accessible to the system with energy E and hence the entropy given by 

S{N,E) = kB\nn{N,E), (3) 

where ks is the Boltzmann constant. Q{N, E) depends on whether the particles are dis- 
tinguishable or not, and whether they are bosons or fermions. Figure [T] shows the single 
particle spectrum. This model is non-trivial because of the presence of degeneracies, for 
example, £(1,0) = e{0, 1) = ^o, and because the energy levels are not simply spaced. This 
makes computing Q{N, E) a non trivial task. Later, to show the versatility of our method, 
we also present results for spin 1 particles and for a system of quantum harmonic oscillators. 
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FIG. 1. The single particle energy spectrum e{p,q) = eo{p'^ + q'^) where p and q are non- negative 
integers. 

B. 'Pen-on-paper' solution 

To start with, it is useful to list the system microstates for a small number of particles. 
For example, for A^ = 3, and E = 20^0, the microstates can easily be listed. Let {pi, qi) 
be the quantum labels for the i*^ particle, then the system microstates {pi,qi',P2,(l2]P3,(l3) 
for indistinguishable spinless particles (bosons) are (2,2;0,0;0,0), (1,2;1,1;0,1), (1,2;1,1;1,0), 
(4,1;1,1;0,1), (4,1;1,1;1,0), (4,0;2,0;0,0), (4,0;1,1;1,1), (0,2;2,0;0,0), (0,2;1,1;1,1), (3,1;3,1;0,0), 
(3,1;3,0;0,1), (3,1;3,0;1,0), (3,1;2,1;2,1), (3,0;3,0;1,1), i.e. there are 14 microstates for indis- 
tinguishable particles. It can be easily inferred that there are 66 microstates for distinguish- 
able particles, and 76 microstates for indistinguishable spin | particles (fermions). Going 
through this exercise leaves one with the distinct impression that it is very cumbersome to 
try to solve this problem using 'pen-on-paper' for a large number of particles. An efficient 
computational algorithm is needed. 
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FIG. 2. Graph of Int versus N, where t is the real computational time in seconds and A^ is the 
number of particles. In this graph the data points are plotted for A^ = 2, 3, 4 and 5 particles. 

C. A brute force computational method 

For distinguishable particles, the simplest brute force method that can be devised involves 
N nested do-loops, each over the list of single particle energy levels. This results in a com- 
putational scheme that scales exponentially with the system size. Still, this is an instructive 
method to apply because it gets the student to appreciate the rapid increase in Q{N, E) for 
moderate A^, and the corresponding exponential increase in computational time. Figure |2] is 
a graph of Int versus A^, where t is the real computational timei^ taken to compute Q{N, E) 
for A^ = 2, 3, 4 and 5 and for E = SOOeo- The graph shows that t oc exp(aA^), where a = 4.61 
is the slope of the curve. We conclude that it will take several million years to consider only 
a moderate number of particles! A more practical algorithm is needed. 

D. A recursive solution for Q{N,E) 



For fixed total energy E, we consider a system of A^ distinguishable particles to be 
composed of a subsystem of (A^ — 1) particles with Q{N ~ 1,E — E ) microstates and a 
subsystem of 1 particle with Q{1,E ) microstates as depicted in Fig. |3l In this context. 




N particles with energy E 



FIG. 3. The system comprising N particles with energy E. 

fi(l, E ) are the degeneracies of the single particle levels with energy E < E. Q{N, E) is 
then simply constructed by summing the number of ways in which the subsystem of (A^ — 1) 
particles can configure itself with energy E — E weighted by the degeneracy of the single 
particle level with energy E , i.e. 

E 

Q{N, E) = Y^ Q{N -1,E-E')xQ{1, E'). (4) 

e'=o 

This is the basis of our recursive algorithm, which can be coded in a very straight forward 
manner. Q{1,E ) can be constructed very simply and, from this, Q{N, E) may be calculated 
recursively for A^ = 2, 3, 4, ■ ■ • . 

A plot of entropy per particle versus energy per particle is presented in Fig. H] for A^ = 3- 
— 10 particles, and compared with the equivalent result calculated within the canonical 
ensemble (see Appendix A) which is, by definition, in the thermodynamic limit. For these 
rather low values of A^, the graphs are not smooth, which makes the temperature ill-defined. 
This is consistent with the fact that temperature for an isolated system only really makes 
sense in the infinite system size limit, a point already alluded to in the introduction in the 
context of Eq. ([1]) . 

Using Eq. (jlj), we are able to compute Q{N,E) for relatively large values of A^. Using 
the computational resources at our disposal, the maximum number of particles that we were 
able to consider is A^ ~ 200. The computational time taken is of the order of minutes. For 
A^ greater than this maximum number, the largest number addressable by the computer- 
is exceeded which results in overflow errors. 

Figure [5] shows a plot of entropy per particle versus energy per particle. We notice that 
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FIG. 4. Entropy per particle versus energy per particle for A^ = 3 — 10 distinguishable particles. 
The result in canonical ensemble (A^ — )• oo) is also included to illustrate where this system is 
relative to the thermodynamic limit (A^ -^ oo). The bulleted point depicts the calculated result of 
il(3, 20eo) = 66 as discussed in subsection II B. 

the result for A^ = 200 is still not converged compared with the result in the infinite system 
size limit. The deviations for the entropy per particle are of order 1% at an energy per 
particle of lO^o- An improved algorithm must be found. 

E. Counting in terms of the smallest real number r that is addressable by the 
computepi^ 



A small improvement in the algorithm can be achieved by counting states in terms of r, 
which is the smallest real number that is addressable by the computer rather than in terms 
of unity which is the normal thing to do. This results in a re-normalized number of states 
accessible to the system Q{N, E) where 



n{l,E) = n{l,E) X r, 



and 



n{N,E) = n{N,E) X r 
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FIG. 5. Entropy per particle versus energy per particle calculated for N = 200, 490 and 2000 
compared with the result calculated in the canonical ensemble {N — )• oo). 

The advantage of counting in terms of r is that a higher value of A^ is achievable before 
Cl{N, E) exceeds the largest number addressable by the computer. Using this method, the 
maximum number of particles that we were able to consider is A^ ~ 500. The computational 
time taken is of the order of minutes. We notice that the result for A^ = 490 in Fig. [5] is 
improved over the value for A^ = 200, but is still not converged compared with the result in 
the infinite system size limit. The deviations for the entropy per particle are of order 0.5% 
at an energy per particle of lO^o- In the next subsection, we describe an algorithm that 
resolves this problem exactly. 

F. A recursive solution for the entropy ^(A^, E) 



The solution to this problem comes from the fact that we are ultimately interested in 
calculating the entropy S'(A^, E) rather than Q{N, E) itself. Rearranging Eq. (|4]) in the 
following manner 

niN, E) = n{N-l, E) fn(l, 0) + ^1!!~^'^~^^ X fi(l, 1) + • • ■ + ^^Z'^'^l X ^(1' E) 



n{N -i,E) 



n{N -i,E) 



(7) 
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FIG. 6. The number of microstates 0,(1, E ) for a single particle for < £" < lOOeo- Its discrete 
structure means that the terms in Eq. @ will be non-zero only when ft{l,E ) is non-zero. 



gives us the following recursive algorithm for the entropy 

E 



S{N,E) = S{N -1,E) 



ks In y. 6xp I 

E'=o ^ 



S{N -l,E-E)-S{N -l,E) 



k 



B 



X n{i,E' 



(8) 



Working directly with ^(A^, E) rather VL{N, E) is much more manageable from a com- 
putational point of view. For the range of A^ and E that we considered in our work, we 
never ran into problems exceeding the maximum number addressable by the computer. The 
computational time taken for N = 2000 particles and E = 2 x W^Eq is of the order of a day. 
We notice that the result for N = 2000 in Fig. [5] is converged compared with the result in 
the infinite system size limit. The deviations for the entropy per particle are less than 0.1% 
at an energy per particle of lO^o- 

G. Microstate complexity versus averaged macrostate properties 

The complexity of this problem is further illustrated by considering the individual terms 
in the summation of Eq. (jl]). We have plotted in Fig. [6l Q[1,E) as a function of E . 
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FIG. 7. A plot of n{N - 1,E - E')/Vl{N - l,E) versus E' (black curve) and 
Q,{N -1,E - E')/n{N -1,E) X n{l,E') versus E' (red curve with circles) where < E' < 
2 X lO^eo and N = 2000 particles. 



This is obviously a discrete curve since Q{1,E ) are the degeneracies of the single particle 
levels with energy E . We have normalized Q{N — 1,E — E ) by Q{N — 1,E) for numer- 
ical convenience and plotted in Fig. [7] for A^ = 2000, n{N -1,E - E')/n{N - 1,E) and 
[n{N -1,E - E')/n{N - 1, E)] X fi(l, E') as a function of E' ior E = 2 x 10%o- Not sur- 
prisingly, the curve of Q{N — 1,E — E ) versus E is smooth and it reflects the convergence 
in the thermodynamic limit of all intensive macrostate variables such as the entropy per 
particle, the chemical potential, the temperature and so on. 



However, Q{N, E) is constructed by summing over terms of the form of the product of 
Q{N — 1,E — E ) with Q{1,E ). These contributions are jagged, distinctly discontinuous and 
rapidly decreasing as a function of E as seen in Fig. [71 The complexity shown in this figure 
is non trivial and shows how the microstate structure of the system fundamentally underpins 
its macrostate properties, even in the thermodynamic limit. No mean- field approach can 
adequately capture this complexity on the microstate scale. 
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H. The heat bath 

It is useful to note that the only contributions to Q{N, E) come from terms for which 
i7(l, E ) are non-zero. So, the subsystem of (A^ — 1) particles can only attain an energy of 
E — E ii Q{1,E ) is non-zero. This is true for any value of A^ and, in particular, this is 
true as A^ — )■ oo. One can therefore be forgiven for thinking about this as a case of the tail 
wagging the dog! 

However, on more careful inspection of Fig. [TJ we note that the main contributions to the 
summation of Eq. (jl]) come from terms close to E = for which the subsystem of (A^ — 1) 
particles has an energy close to E. So, all but the low energy single particle states are 
suppressed. In the large A^ limit, one may view the system of A^ particles as being composed 
of one particle in thermal equilibrium with the (A^ — 1) particle subsystem; the latter may 
therefore be viewed as a heat bath. This picture enables us to arrive at the Boltzmann 
distribution that cements the relation with the canonical ensemble. These observations 
have important implications for the determination of the temperature of the system, which 
we consider more completely next. 



I. Temperature 

A simple way to determine the temperature of the system as a function of energy per 
particle is to apply Eq. ([T]) to the curve of entropy per particle versus energy per particle 
in the large A^ limit. To demonstrate the versatility of our methods we have considered a 
different model now: We have plotted the temperature of a system of distinguishable spin 1 
particles as a function of A^ and at energy per particle fixed at 0.5eo in Fig- El In the limit 
of large A^ the numerical result converges to the analytical result as it should since these 
results are derived in the limit of large A^. The percentage difference between the result in 
the canonical ensemble (A^ — )■ oo) and the numerical result for A^ = 5000 spin 1 particles is 
0.03%. 

A more interesting way to extract the temperature that cements the equivalence of the 
microcanonical ensemble with the canonical ensemble is now presented. Following Eq. (|8]), 
and for A^ = 2000, which we have already demonstrated is sufficiently close to the thermo- 
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FIG. 8. A plot of the inverse temperature as a function of the number of spin 1 particles, A'^, 
at fixed energy per particle equal to O.Seo- The canonical result {N — )• oo limit) for the inverse 
temperature at the same energy per particle is also plotted for comparison. Analytical expressions 
for thermodynamic quantities within the microcanonical ensemble such as the temperature of the 
spin 1 system are usually obtained using the method of Lagrange multipliers.—*^ 

dynamic limit, we have plotted in Fig. [9l 



S{N -1,E-E') - S{N - 1, E) 



k 



B 



(9) 



as a function of i? for -E = 2 x lO^eo- 

The results are negative, as they should be, and linear for low values of E . The curve 
deviates downwards, i.e. becomes more negative, for higher values of E which corresponds to 
the super-suppression of the higher energy single particle states. This super linear behaviour 
for higher values of E we found to be very intriguing and we are not aware that this has 
been noted before. For the low energy single particle states, we may adequately model this 
curve by 



A = -PE 
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FIG. 9. A plot of A = S{N-i,E~E^)^~s{N-i,E) ^g ^ function of E' iov E = 2 x lO^eo and for 
< E < 2 X lO^eo- The curve is linear for low vales of E . The broken curve in this plot is the 
linear region extrapolated to i? = 2 x lO'^eo to show the downward deviation of the curve of A for 
higher values of E , which corresponds to the super-suppression of the higher energy single particle 
states. The linear region of the curve of A versus E has a slope equal to — 0.0638e,^ , which we 
identify as — /3 according to Eq. (jlOp . 



so that Eq. ([H]) may now be re-written as 

oo 

S{N, E) = S{N ~l,E) + kB In ^ exp(-/3 E') x f](l, E'] 



(11) 



E =0 



This works fine for the higher energy states as well which are suppressed even further than 
what the linear relation in Eq. (TTOj) implies. The upper limit in the summation in Eq. ( TTTj) 
has therefore, accordingly, been replaced with oo with no change to the final result. /3, of 
course, is a function of E and emerges here simply from the linear fit to the low energy 
data in Fig. O This establishes the equivalence with the canonical ensemble for which /3 is 
identified with the inverse temperature 



P 



1 

k^ 



(12) 
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and z{T) with the single particle partition function 

oo 

J2 exp{-f3E') X n{l,E'). (13) 



oo 
Z = 

E'=0 



J. Chemical potential 

Since Eq. (jlj) corresponds to us systematically increasing the system size by one particle 
at a time, a very natural thermodynamic quantity to consider is the chemical potential which 
is defined as 

Since it is very difficult to numerically keep the entropy fixed (the preceding sections will 
attest to this!), it is more useful for us to consider the quantity 

instead. Because 6N = 1, we may write 

-^ = S{N,E)-S{N-1,E), (16) 

which can be easily evaluated using Eq. ([H]). 

Plotting — ^ as a function of A^ will give an indication of the convergence of — ^ in 
the thermodynamic limit. However, there is a subtle point that is worth noting. The 
differentiation process in Eq. (fTSjl requires that the total energy E must be kept constant. 
However, if we keep E constant in Eq. flT6l) . then E/N ^' as A^ — ?■ oo. In this case, the 
entropy per particle — )■ as A^ — )■ oo, and Eq. (IT^ simply converges to zero. This is correct 
but not very useful! It is more instructive to keep the energy per particle fixed in Eq. fITBl) . 
and in so doing the correct convergence properties of — ^ will be achieved. We leave this as 
an exercise for the reader to investigate. 

A second subtle point related to the chemical potential as defined in Eq. flT6|) is that 
the left hand side of this equation is manifestly intensive, whereas the right hand side is 
manifestly extensive. The reader is encouraged to resolve this dilemma. The results can be 
investigated numerically according to the discussion in the previous paragraph. 

Another useful way to approach an investigation of the chemical potential is to consider 
the following scaling argument. Let s{e) be the entropy per particle in the thermodynamic 
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FIG. 10. A plot {—[ijT) as a function of the number of harmonic oscillators, N ^ at fixed energy 
per particle equal to lOeo- The canonical result (A^ — )• oo limit) for the chemical potential at the 
same energy per particle is also plotted for comparison. The analytical expression for the inverse 
temperature as a function of the energy per particle, e, is — /^/T = ks ln(e + 1). 






limit, where e is the energy per particle. Then 



S{N,E) = Ns 



Applying Eq. (IT311 to the above expression gives 



/i , , ds(e) 



(17) 



(18) 



Once again, to demonstrate the versatility of our method, we have plotted this expression 
in Fig. [To] for a system of quantum harmonic oscillators as a function of A^ to test the 
convergence of — ^ in the thermodynamic limit. This is done at constant energy per particle 
e equal to lOeo rather than constant total energy E. The numerical result for the chemical 
potential calculated using Eq. ( ITSl) converges to the analytical result in the limit of large A^. 
The percentage difference between the analytical result (in the limit of A^ — ^ oo ) and the 
numerical result for A^ = 2000 harmonic oscillators is 0.1%. 
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III. FINAL REMARKS 

We have developed a hierarchy of methods that enable the computation of the entropy 
for a system of independent and distinguishable particles. Our methods are very general and 
applicable to a very wide class of problems. Our algorithm expressed in Eq. ([H]) is robust and 
may be studied in the large A^ limit with a reasonable computational effort. To within limits 
set by the computer, A^ can be made arbitrarily large. This method enables one to monitor 
the approach to the thermodynamic limit as a function of A^ for various thermodynamic 
quantities such as the entropy per particle, the temperature, the chemical potential, and so 
on. 

We have given some exercises for the reader to pursue this subject more completely. It 
is helpful to first apply our expressions to an analytically solvable model such as a system 
of independent quantum harmonic oscillators or a system of independent spin 1 particles. 
Some of the arguments, for example those referred to the chemical potential in Section II J, 
may first be tested against known analytical results. 

Considering indistinguishable spinless particles (bosons) or indistinguishable spin | par- 
ticles (fermions) would be a natural extension to our methods. 

Appendix A: The canonical ensemble 

In the text, results for various properties in the microcanonical ensemble for finite N are 
compared with equivalent results in the canonical ensemble which are, by definition, in the 
thermodynamic limit. The relevant expressions for the canonical ensemble are given below. 

The single-particle partition function z{X) for the model described in Eq. ([2]) is given by 

oo oo 

z{X) = ^exp(-Xp2) X 5^exp(-Xg^), (Al) 

p=0 q=0 

where p and q are non-negative integers, and where 

is dimensionless and is a measure of the inverse temperature. z{X) as a function of X can 
be computed numerically in a very straight-forward manner. 

From this, the Helmholtz free energy per particle f{X) in units of Eq is given by 

f{X) = -^lnz(X), (A3) 
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and the average total energy per particle e{X) in units of eo is given by 

e(X) = -^ln^(X). (A4) 

We calculate the average entropy per particle in units of ks from 

six) =X{e{X) -fix)). (A5) 

This enables us to plot, for example, the entropy per particle versus the energy per particle 
(at the same value of X). 

Another useful result for our analysis is the expression for the chemical potential in units 
of £0 which is given by 

/i = fix). (A6) 
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